arXiv:1504.07508v3 [cond-mat.mtrl-sci] 23 Sep 2015 


Lattice dynamics and vibrational spectra of the orthorhombic, tetragonal and cubic 

phases of methylammonium lead iodide 

Federico Brivio,^ Jarvist M. Frost/ Jonathan M. Skelton/ Adam J. Jackson/ Oliver J. Weber/ Mark 
T. Weller/ Alejandro R. Goni/ Aurelien M. A. Leguy/ Piers R. F. Barnes/ and Aron Walsh^A.j^ 

^ Centre for Sustainable Chemical Technologies and Department of Chemistry, 

University of Bath, Claverton Down, Bath BA2 7AY, UK 
^ICREA, Passeig Lluts Companys 23, E-08010 Barcelona, 

Spain; Institut de Ciencia de Materials de Barcelona (ICMAB-CSIC), Campus UAB, E-08193 Bellaterra, Spam 
^Department of Physics, Imperial College London, SW7 2AZ, UK 
‘^Global E3 Institute and Department of Materials Science and Engineering, Yonsei University, Seoul 120-749, Korea 

(Dated: September 24, 2015) 

The hybrid halide perovskite CHjNHgPblg exhibits a complex structural behaviour, with succes¬ 
sive transitions between orthorhombic, tetragonal and cubic polymorphs at ca. 165 K and 327 K. 

Herein we report first-principles lattice dynamics (phonon spectrum) for each phase of CHgNHgPbIg. 

The equilibrium structures compare well to solutions of temperature-dependent powder neutron 
diffraction. By following the normal modes we calculate infrared and Raman intensities of the 
vibrations, and compare them to the measurement of a single crystal where the Raman laser is 
controlled to avoid degradation of the sample. Despite a clear separation in energy between low 
frequency modes associated with the inorganic (Pblg),^ network and high-frequency modes of the 
organic CHjNH/ cation, significant coupling between them is found, which emphasises the inter¬ 
play between molecular orientation and the corner-sharing octahedral networks in the structural 
transformations. Soft modes are found at the boundary of the Brillouin zone of the cubic phase, 
consistent with displacive instabilities and anharmonicity involving tilting of the Pbig octahedra 
around room temperature. 

PACS numbers: 63.20.D,63.20.Ry,78.30.-j 


I. INTRODUCTION 

Materials that adopt the perovskite crystal structure 
are known for their complex structural landscapes, with 
a large number of accessible polymorphs depending on 
the temperature, pressure, and/or applied electric held. 
For ternary ABXg perovskites, the A site cation is at the 
centre of a cube formed of corner sharing BXg octahedra. 
Displacement of the A cation is usually associated with 
a ferroelectric (Brillouin zone centre) instability, while 
tilting of the BXg octahedral network is usually linked to 
antiferroelectric (Brillouin zone boundary) transitions.^^ 

Hybrid organic-inorganic perovskites are formed when 
one of the elemental perovskite buildi ng blocks is re¬ 
placed by a molecular anion or cation.l^ There exists 
a large family of such compounds, including the widely 
studied formate perovskites, which contain both molecu¬ 
lar anions and cations.^^^t^ Hybrid halide perovskites are 
of current intense research effort, due to their high effi¬ 
ciency photovoltaic action.I^lti^ 

Methylammonium lead iodide (MAPbIg, where MA 
represents the CHgNHj/ cation), was hrst reported by 
Weber in 1978.1^^ It is the most relevant hybrid halide 
perovskite for photovoltaic application. The transition 
from orthorhombic to tetragonal to cubic perovskite 
structures as a function of temperature has been stud¬ 
ied by techniques including calorimetry and infrared 
spectroscopjff^, single-crystal X-ray diffractioiP^, and di¬ 
electric spectroscopy.!^ Recently analysis of powder neu¬ 
tron diffraction (PND) measurements has provided more 


quantitative insights into the temperature dependent be¬ 
haviour of the MA cation within the anionic (Pblg),^ 
network.!^ There is now direct evidence for the degree 
of order of the MA cation in the different phases, and 
the average lattice parameters (and thus extent of octa¬ 
hedral tilting) as a function of temperature through the 
first and second order phase transitions. Quasi-elastic 
neutron scattering has provided further insights into the 
rotational dynamics of the MA cation with a room tem¬ 
perature residence time of ~ 14 ps,!^ while time-resolved 
vibrational spectroscopy identified fast librations (300 fs) 
and slow (3 ps) rotations of the molecule.!^ 

In this study, we calculate the phonon dispersion in 
each phase of MAPbIg within the harmonic approxima¬ 
tion, computing the force constants with density func¬ 
tional theory (DFT). We use the PBEsol functional, 
which is a generalised-gradient-approximation (GGA) to 
the exchange-correlation functional, numerically evalu¬ 
ated with Perdew’s method, adjusted to give more ac¬ 
curate lattice constants and forces for solids.!^ The lat¬ 
tice dynamic calculations allow the atomic origin of each 
phonon mode to be identified. Changes in lattice po¬ 
larisation and polarisability for each eigenvector provide 
the infrared and Raman activity of each mode. Spec¬ 
tral features related to the inorganic and organic com¬ 
ponents (from 0 to 3000 cm“^) are well reproduced in 
comparison to the Raman spectra of a single crystal of 
MAPbIg. Overlap is found between the vibrations of the 
CHgNH/ and PbIg components up to 130 cm“^, with the 
modes from 300 to 3000 cm“^ being associated with pure 
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molecular vibrations. The phonon dispersion has impli¬ 
cations for developing quantitative models for the gen¬ 
eration, transport and recombination of photo-generated 
electrons and holes in hybrid perovskite solar cells. 

A. Structure Models 

The normal modes of a system are defined for an 
equilibrium configuration. Calculating the vibrations 
for a non-equilibrium structure will result in imaginary 
frequencies upon diagonalising the dynamical matrix. 
Therefore we have generated well optimised structures 
of MAPbl 3 . One challenge in calculating the phonons 
of hybrid perovskites is the soft nature and complicated 
potential energy landscape of some of the restoring po¬ 
tentials, particularly those involving the organic cation. 

The models for the crystal structures used in this study 
are discussed in detail below and a comparison with the 
measured diffraction patterns are provided as Supple¬ 
mentary Material.l^ 


of quasi-thermal equilibrium. It is possible that differ¬ 
ent ordering might be frozen into the low temperature 
phase by epitaxy or application of external force or elec¬ 
tric fields. 



1. Orthorhombic Phase 




The orthorhombic perovskite structure is the low tem¬ 
perature ground state of MAP bl.^ and maintains its sta¬ 
bility up to ca. 165 A comparison of the en¬ 

thalpy from DPT calculations confirms this ordering in 
stability. The difference in enthalpy is small, just 2 meV 
per MAPbIg unit compared to the most stable tetrago¬ 
nal phase, yet 90 meV compared to the high-temperature 
cubic phase. 

Initial dif fractio n pattern solutions assigned the Pna2i 
space group.li^I^ Recent analysis of higher quality pow¬ 
der neutron diffraction data reassigns it to Pnma (a D 2 h 
point group).The structure is a v^a x \f2a x 2a su¬ 
percell expansion of the simple cubic perovskite lattice, 
i.e. following the lattice transformation matrix 


In the Pna2i phase, the Pbig octahedra are distorted 
and tilt as a^b~b~ in Glazer notatioiP with respect to 
the orientation of the conventional cubic cell. In this 
low-temperature phase, the four molecular cations in the 
unit cell are static on the diagonals of the ab planes point¬ 
ing towards the undistorted facets of the cuboctahedral 
cavity. Correspondingly, molecules belonging to differ¬ 
ent planes are anti-aligned with a head-tail motif. Such 
an antiferroelectric alignment is expected from consider¬ 
ation of the molecular dipole-dipole interaction.^^ 

In the low temperature orthorhombic phase the 
CH 3 NH 3 sublattice is fully ordered (a low entropy state). 
The ordering may be sensitive to the material prepara¬ 
tion and / or cooling rate into this phase, i.e. the degree 


FIG. 1. (Color online) The crystal structures of the 
(a) orthorhombic, (b) tetragonal and (c) cubic phases of 
CHgNHgPbIg. Thc upper and lower panels are oriented 
through < 100 > and < 001 >, respectively. Lattice param¬ 
eters and coordinates obtained from powder neutron diffrac¬ 
tion were optimised using density functional theory (PBEsol). 
The PbIg octahedra are shaded grey. All structures are avail¬ 
able in an on-line repository.^ 


2. Tetragonal Phase 

At 165 K MAPbIg goes through a first-order phase 
transition from the orthorhombic to the tetragonal space 
group lA/incm {D^h point group), which continuously 
undergoes a second- order p hase transition to the cubic 
phase by ca. 327 k | 21 | 26 | 27 | orthorhombic 

phase, this can be considered a ^/2a x ^/2a x 2 a expansion 
of the cubic perovskite unit cell. 

The molecular cations are no longer in a fixed posi¬ 
tion as in the orthorhombic phase. The molecules are 
diso rdered between two non-equivalent positions in each 
cage 122130] tetragonal distortion parameter in the cu¬ 
bic basis is greater than unity ~ 1.01 at 300 K), cor¬ 
responding to an elongation of the Pbig octahedra along 
the c axis. The associated octahedral tilting pattern is 
aPaPc~ in Glazer notation. 

Atomistic simulations within periodic boundary condi¬ 
tions require an ordered configuration. The solved crys¬ 
tal structure shows that there are several possible con¬ 
figurations for the organic cations within the tetragonal 
unit cell. These configurations have similar enthalpies 
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Phase a (A) b (A) c (A) d(Pb-I) (A) Z Cut-off (eV) A:-points Forces (meV/A) AH (meV) 


Orthorhombic 










DFT/PBEsol 

9.04 

12.66 

8.35 

3.18 

4 

700 

5x4x5 

1 

0 

PND (100 K) 

8.87 

12.63 

8.58 

3.19 






Tetragonal 










DFT/PBEsol 

8.70 

8.72 

12.83 

3.19 

4 

800 

5x5x3 

1 

2 

PND (180 K) 

8.81 

8.81 

12.71 

3.17 






Cubic 










DFT/PBEsol 

6.29 

6.23 

6.37 

3.17 

1 

700 

6x6x6 

1 

90 

PND (352 K) 

6.32 

6.32 

6.32 

3.16 







TABLE I. Equilibrium cell parameters from DFT/PBEsol energy minimisation, including the converged plane-wave cut-off, 
fc-point mesh and force threshold. Z represents the number of formula units of CHgNHgPbIg per cell. The calculated difference 
in enthalpy (AH) of each phase is given with respect to the ground-state orthorhombic configuration and per CHgNHgPbIg unit. 
Shown for comparison are the cell parameters and average Pb-I inter-atomic separations (d) from powder neutron diffraction 
(PND).>^ 


within DFTjl^ which is consistent with the observed dis¬ 
order. We choose to use the most energetically stable 
structure, which is also consistent with a previous DFT 
investigation.!^ 

In the model of the tetragonal structure, the MA 
cations are aligned as in the orthorhombic phase, to¬ 
wards the face of the perovskite cage, i.e. < 100 > in 
the cubic basis. The MA in different (001) planes are ap¬ 
proximately orthogonal to one other. The orientional dy¬ 
namics of the methylammonium ions, which exists above 
165 K, is not taken into account in this equilibrium con¬ 
figuration study. 


3. Cubic Phase 

With increasing temperature the tetragonal lattice pa¬ 
rameters become more isotropic (i.e. ^ moves closer 

to 1), and the molecular disorder increases, to the point 
where a transition to a cubic phase occurs around 327 
K. The transition can be seen clearly from changes in 
the heat capacity,!^ as well as in temperature dependent 
neutron diffractiorP^!. 

The cubic space group PmZm {Oh symmetry) has been 
assigned to this high-temperature phase. Although the 
methylammonium ions posses symmetry, the orien¬ 
tational disorder gives rise to the effective higher symme¬ 
try on average. The local structure will necessarily have 
a lower symmetry. Indeed, for the bromide and chloride 
analogues of MAPbIg, pair-distribution function analy¬ 
sis of X-ray scattering data indicates a local structure 
with significant distortion of the lead halide framework 
at room temperature.!^ 

We previously considered alignment of the molecules 
along the principal < 100 > (face), < 110 > (edge), and 

< 111 > (diagonal) directions of the cubic unit cell, and 
showed that they are of similar DFT enthalpy, with a 
small barrier for rotation.!^ Further ab-initio molecular 
dynamics showed an average preference for the < 100 > 
facial configuration at 300 Therefore we chose the 

< 100 > configuration as our reference structure for the 


lattice vibrations. 

Representations of the crystal structure of each phase 
are shown in Figure the equilibrium structure param¬ 
eters are listed in Table and the structures themselves 
are available in an on-line repository.!^ 


II. METHODS 
A. Computational 

The total energy and atomic forces were computed 
from first-principles within den sity f unctional theory as 
implemented in the code VASP.!^3^ Noise in the lattice 
vibrations was minimised by rigorous convergence of total 
Kohn-Sham energy with respect to the basis set (kinetic 
energy cut-off for plane waves) and sampling of reciprocal 
space (density of the fc-point mesh). The final computa¬ 
tional set-up is summarised in Table |l] 

We performed complete optimisation of the cell vol¬ 
ume, shape and atomic positions, with the PBEsoP^ 
semi-local exchange-correlation functional. The scalar- 
relativistic projector-augmented wave methocP^was em¬ 
ployed, with a pseudo-potential treating the Pb 5d or¬ 
bitals as valence. Spin-orbit coupling was not considered 
as it mainly affects the Pb 6p conduction band, which 
does not influence the interatomic interactions at equilib¬ 
rium. All atomic forces were reduced to below a threshold 
of 1 meV /A. Due to the presence of the organic cations, 
which breaks the ideal lattice symmetry, deviations in the 
expected parameters can occur, e.g. in the high temper¬ 
ature pseudo-cubic phase, the three equilibrium lattice 
parameters are not equal. The equilibrium structure pa¬ 
rameters (at 0 K and excluding zero-point contributions) 
are reported in Table 

The normal modes are calculated within the har¬ 
monic approximation, using the PhonopjEMSI package 
to construct and evaluate the dynamical matrix com¬ 
posed of DFT force constants. Both the finite displace¬ 
ment method (FDM or supercell approach)^ and density 
functional perturbation theory (DFPT)!^ approaches to 
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FIG. 2. (Color online) (a-c) Projected phonon density of states (PDOS) for the three phases of CHgNHgPbIg as calcnlated 
from DFT/PBEsol, generated by convolution with a 32 cm“^ Lorentzian. Note that no imaginary (negative) modes are found 
in the orthorhombic or tetragonal phases, (d-f) Simulated infrared (IR) spectra, (e) The measured IR spectrum at 300 K 
from Ref. [35] is coplotted in mustard (light grey), (g-i) Simulated Raman spectrum, (g) The measured Raman spectrum of 
a single crystal in the orthorhombic phase at 100 K is coplotted in orange (light grey). Simulated spectra were broadened 
by convolution of a 2 cm“^ Lorentzian. Insets in the simulated spectra are the same data on a logarithmic scale to show the 
structure in the low intensity modes. 


construct the force constants were tested. The results 
of both approaches produced similar vibrational spectra, 
with a variance in the mode energies of 6 cm“^. 

Within a primitive cell of JV atoms there are 6N pos¬ 
sible displacements (iLx, iLy, iLz), which can be reduced 
by the crystal symmetry. For the orthorhomic phase, the 
288 possible displacements are reduced to 41, while the 
tetragonal and cubic phases required 288 and 72 displace¬ 
ments, respectively. The phonon dispersion (for q-points 


away from the Brillouin zone centre, the F point) in the 
cubic phase was probed in a 2 x 2 x 2 supercell. Due to 
computational expense, we do not calculate this for the 
other (larger unit cell) phases. For the phonon density 
of states, Brillouin-zone integrations were performed with 
36 X 36 X 36 (orthorhombic and tetragonal) and 48 x 48 x 48 
(cubic) F-centered Monkhorst-Pack g-meshes. 

Once the normal eigenmodes and eigenvalues are calcu¬ 
lated, it is possible to model their associated Raman and 
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infrared (IR) activity by mode following. The two spec¬ 
troscopic techniques probe different physical responses of 
the material: the change in polarisation for IR and the 
change in polarisability for Raman. The IR spectra are 
simulated with the analytic formula of Gianozzi & Baroni 
(using the Born effective charge tensorjl^^. Prediction 
of the Raman spectra required computing the change in 
macroscopic dielectric tensor with respect to each nor¬ 
mal mode of the system, a significant DFT calculation in 
terms of computational expense.!^ 


B. Experimental 

Methylammonium lead iodide single crystals were 
grown according to the method of Poglitsch and Weber.l^ 
12.5 g of lead acetate trihydrate (Pb(CH3C02)2 • 3 H 2 O, 
Sigma) was dissolved in 10 mL hydroiodic acid (Hlaq, 57 
wt%, Sigma) in a 50 mL round bottom flask and heated 
to 100°C in an oil bath. Separately, 0.597 g of CH 3 NH 2 
(aq, 40 wt%, Sigma) was added dropwise to a further 2 
mL of Hlaq kept at 0°C in an ice bath under stirring. The 
methylammonium iodide solution was then added to the 
lead acetate solution and the mixture was cooled over five 
days to a temperature of 46°C, resulting in the formation 
of black crystals with largest face length around 8 mm. 
The content of the flasks was subsequently filtered and 
dried for 12 hours at 100°C. 

Raman spectra were collected in backscattering geom¬ 
etry with a high resolution LabRam HR800 spectrom¬ 
eter using a grating with 600 lines per millimetre and 
equipped with a liquid-nitrogen-cooled charge coupled 
device (CCD) detector. The 785 nm line of a diode- 
pumped solid state laser was used as excitation beam 
and focused onto the sample using a long distance 20 x 
microscope objective. Raman measurements were carried 
out at 100 K using a gas-flow-type cryostat with optical 
access that fits under the microscope of the Raman setup. 
The high spectral resolution and stray-light rejection of 
the LabRam spectrometer, particularly in combination 
with the 785 nm line, allowed us to measure the Raman 
spectrum of MAPI at very low Raman shifts down to 
20 to 30 cm“^. This way and for the first time, we were 
able to spectrally resolve several low-frequency modes as¬ 
sociated to vibrations of the inorganic cage of the hybrid 
perovskite. 

Heating by laser light directly absorbed by 
CH 3 NH 3 Pbl 3 has been shown to lead to rapid 
degradation of the material resulting in Pbl 2 Ra¬ 
man signatures.!^ Since 785 nm light is only weakly 
absorbed the heating effect of the laser was low enough 
to ensure the crystal structure was preserved. The power 
density incident on the sample was kept at 80 W/cm^. 
At such power level, it was checked that no appreciable 
spectral changes in peak width and/or position occurred, 
yet maintaining a good signal-to-noise ratio. Further, 
samples were kept under vacuum inside the cryostat 
during the measurements. 


III. RESULTS 
A. Harmonic Phonons 

The full phonon density of states (DOS) is shown for 
the three phases of CH 3 NH 3 Pbl 3 in Figure]^ Also plot¬ 
ted is the partial DOS, where assignment to CH 3 NHj^ or 
Pbljj is performed based on the atomic contribution to 
each eigenvector. An animation of all 36 eigenmodes of 
the cubic phase is provided as Supplementary MateriaP^. 

Qualitatively, each MAPbl 3 phase shows similar vibra¬ 
tional properties with three energetic regions of phonons: 
(i) a low-frequency band from 0-150 cm“^; (ii) a 
mid-frequency band from 280-1600 cm“^; (iii) a high- 
frequency band from 2900-3300 cm“^. These r anges are 
consistent with previous computational reports.!^^*^ 

Due to the large difference in atomic mass of the or¬ 
ganic and inorganic components, and to the difference in 
bonding between the inorganic cage and the covalently 
bonded molecule, we anticipated that the low frequency 
modes would comprise entirely of motion by the Pbig oc- 
tahedra, while the high frequency modes will involve the 
CH 3 NH;/ cation. This is qualitatively the case, but there 
is also significant coupling between the two. Taking the 
example of the cubic phase with 36 modes: the highest- 
energy 18 modes (forming bands ii and iii) correspond to 
molecular vibrations, i.e. the 3A-6 modes of the methy¬ 
lammonium ion. For an isolated non-linear molecule, the 
6 translational and rotational degrees of freedom do not 
contribute to the pure vibrational spectrum, but this is 
not the case for a molecule inside a cuboctahedral cavity. 

The 6 additional molecular modes are strongly cou¬ 
pled to the 9 modes (i.e. 3A^-3 for PbIg) associated with 
stretching of the Pb-I bonds and breathing of the Pbig 
octahedra, which results in the spectral overlap observed 
in the partial DOS of band i. Particularly striking is the 
low-frequency pivoting motion associated with the libra- 
tion of the molecular dipole, coupled with a breathing of 
the octahedral framework (e.g. modes 10 and 15 in the 
SI). The final three zero-frequency modes correspond to 
acoustic translations of the lattice. 


B. Vibrational Spectra 

For MAPbIg it is not possible for the relaxed (non- 
idealised) crystal structures to assign the spectral activ¬ 
ity directly with group theoretic irreducible representa¬ 
tion analysis of the phonon modes. The molecule breaks 
the average crystal symmetry. This symmetry lowering 
allows for simultaneous Raman and IR activity even in 
the pseudo-cubic phase. The predicted spectra (the F 
point phonon modes weighted by the computed spectral 
intensity, convolved with a Lorentizian for experimental 
comparison) are reported in Figure]^ for each phase. 

Raman and IR activity is observed across each of the 
three phonon bands previously discussed. A notable ex¬ 
ception is the lowest-energy purely molecular vibration 
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FIG. 3. (Color online) Comparison of the calcnlated and measured (100 K) Raman spectrnm for the orthorhombic perovskite 
phase of CH 3 NH 3 Pbl 3 in the range 0-1800 cm“^, which includes phonon modes centred on the inorganic cage (lowest frequency), 
organic cage (highest frequency) and coupled-modes (intermediate frequency). The simulated spectrum was broadened by 
convolving a 2 cm“'^ Lorentzian (to match experimental broadening of the high resolution Raman spectrometer), while the 
underlying phonon density of states is shown for comparison with a broadening of 0.5 THz. 


near 300 cm which is neither Raman nor IR active. 

To understand the effect of embedding the methylam- 
monium in MAPbl 3 , we calculate the normal mode vi¬ 
brations of an isolated methylammonium ion in vacuum 
with a similar density functional theory method to our 
periodic calculations (PBE with an atom centered aug¬ 
mented cc-pVQZ basis set). Thereby we calculate the 18 
molecular modes, directly accessing their symmetries and 
nature. The C 3 „ symmetry of the molecule separates the 
vibrational bands into one A symmetric mode and blue- 
shifted two-fold degenerate asymmetric E modes. The 
six bands we find are in ascending energy: twist around 
the C-N axis (282, 886 cm“^); vibration along the C- 
N axis (923, 1239 cm“^); bending of the C-H bonds 
(1418, 1451 cm“^); bending of the N-H bonds (1478, 
1622 cm“^); stretching of the C-H bonds (3018, 3119 
cm“^); stretching of the N-H bonds (3321, 3395 cm“^). 

As the cation charge density is centred towards the 
N, the motion of the protons associated with N have 
the strongest affect on the dipole moment and there¬ 
fore strongest IR activity. Due to the stronger bonds, 
their frequencies are consistently blue shifted relative to 
the C end. Owing to the molecular dipole moment (2.2 
D,[^ which is rotation and position invariant, and corre¬ 
sponds to a polarisation contribution of ~ 3 /iC/cm^), 
the two high-frequency asymmetric stretching modes of 
NHj^ (band iii) results in the strongest absolute IR in¬ 
tensity. The hydrogen stretching modes (band iii) are re¬ 
sponsible for significant Raman activity. The only mode 
involving C or N motion is the weakly IR and Raman 
active vibration at 923 cm“^(vacuum), 1007 cm“^(cubic 
perovskite). 


The rotation of the CH 3 against the NHg unit, while 
being strongly populated in molecular dynamic sim¬ 
ulations, and which forms the main source of quasi¬ 
inelastic neutron scattering, is entirely IR and Raman 
inactive in vacuum. This is the mode responsible 
for the 282 cm“^(vacuum), 318 cm“^(cubic) , 300-310 
cm“^(tetragonal), 370-372 cm“^(orthorhombic) vibra¬ 
tion. Progressive confinement of MA from vacuum to 
the orthorhombic phase blue-shifts the energy of the vi¬ 
bration. 

In the solid state, the degeneracies in the molecular 
modes are typically split by local environment effects, 
peaks are both redshifted and blue-shifted, and the IR 
and Raman activity varies. As such, it is evident that 
analysis of the Raman and IR spectra in the experimen¬ 
tally easily accessible molecular frequency range can en¬ 
able statements to be made about the local structure and 
configuration of the hybrid perovskite. Our data is col¬ 
lected for a particular representation of the cubic and 
tetragonal phase; in reality the location of the MA in 
these phases will be disordered. As such, detailed com¬ 
parison of theory to experiment will require sampling the 
thermodynamic ensemble of structures. 

To our knowledge, no Raman spectra for the three 
phases have previously been reported across the full fre¬ 
quency range. The spectrum up to 450 cm“^ was re¬ 
ported in Ref. Reliable measurements are a challenge 
due to chemical instability of the material. MAPbIg is 
strongly affected by environmen tal co nditions, such as 
the presence of ambient moistur^^ZIlIl Isolated in vac¬ 
uum the material can still decompose and bleach due 
to heating, including by that imposed by the (typically 
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high) Raman laser fluence!^ Such degradation leads to 
the formation of Pbl2, which overlaps in Raman spectra 
with MAPbIg, and so easily leads to misinterpretation. 

The Raman spectra of a high-quality single crystal of 
MAPbIg is shown in Figure and compared in detail 
with the calculations in Figure Across the full spec¬ 
tral range, the agreement between the predicted and mea¬ 
sured spectra is good, with the response across bands i, ii 
and iii well reproduced. On closer inspection in Figure]^ 
there are noticeable shifts in peak positions, which can be 
attributed to three potential sources of error: (a) the har¬ 
monic approximation (anharmonic renormalisation may 
be large); (b) the limits of the exchange-correlation treat¬ 
ment (non-local interactions may be important); (c) the 
assumption of a fully-ordered structure (local inhomogen- 
ity may be prevalent). There is also a notable case of a 
missing peak around 150 cm“^, which we can tentatively 
attribute to a lifetime broadening effect. The same level 
of theory applied to the lead based semiconductors PbS 
and PbTe results in quantitative agreement with mea¬ 
sured phonon frequencies and dispersion,!^ which high¬ 
lights the complex nature of CHgNHgPbIg. 

We have also included in Figure a room tempera¬ 
ture IR spectrum reported by Glaser et al.!^, which again 
shows excellent agreement across the spectrum. A num¬ 
ber of very weak absorption peaks below 3000 cm“^ are 
evident, which could be related to molecular disorder 
and/or partial decomposition. The temperature resolved 
(between 140 and 299 K) IR spectra for tetragonal and 
orthorhombic phases have been previously reported by 
YamamurcpSl. We reproduce the position and the inten¬ 
sity of the peak observed at 900 cm“^, reliably assigning 
it to the (also Raman active) C-N bond stretch. High 
quality IR measurements, in particular looking at the 
very low energy transitions, would provide considerable 
information on the nature of the domains and local struc¬ 
ture in a MAPbIg film. 


C. Anharmonic Effects 

The lattice dynamic simulations discussed above were 
performed within the harmonic approximation. All 
eigenmodes at the centre of the Brillouin zone were real 
(positive frequencies) for each phase, i.e. the structures 
are locally stable. 

The phonon dispersion across the first Brillouin zone 
is shown for the cubic perovskite structure in Figure 
Here imaginary (negative frequency or ‘soft’) modes are 
found at the zone boundaries. Such instabilities are a 
common feature of the perovskite structure, and repre¬ 
sent antiferroelectric distortions linked to to rotations 
and tiling of the octahedra in neighbouring unit cells.!^ 
The soft modes are centred around the R and M points, 
which correspond to the < 111 > and < 110 > direc¬ 
tions in the cubic lattice. This behaviour is similar to 
the inorganic perovskite CsPbClg, where neutron scat¬ 
tering was used to probe condensation of these modes. 





FIG. 4. (Color online) Vibrational properties of the cnbic 
phase of MAPbIg: (a) Phonon dispersion of the low freqnency 
(inorganic cage) modes within the harmonic approximation. 
Negative frequency (imaginary or ‘soft’) modes are found 
at the Brillouin zone boundary points M (^(|i|,0)) and 
R (^(|, I)), (b) Volumetric thermal expansion within the 

quasi-harmonic approximation, (c) Average Griineisen pa¬ 
rameter within the quasi-harmonic approximation. 


which leads to successive transitions from the cubic to 
tetragonal to orthorhombic phases.!^ The effect of these 
modes in MAPbIg, and the associated high levels of an- 
harmonicity at room temperature can be observed di¬ 
rectly in molecular dynamics simulations, where tempo¬ 
ral rotations of the CHgNHj^ ions and distortion s of the 
Pbig octahedra have been found in several studies 

An approach to including the effects of temperature 
(thermal expansion) and first-order anharmonicity in lat¬ 
tice dynamic c alcula tions is the quasi-harmonic approxi¬ 
mation (QHA) .IS21S4I The computational cost is one order 
of magnitude higher than the harmonic approximation 
and thus was considered for the cubic phase only. The 
volumetric thermal expansion coefficient extracted from 
the PND data at 300 K is 1.32 x I0“‘*/K, which compares 
very well to the value of 1.25 x 10“^/K computed within 
the QHA. The predicted thermal expansion for MAPbIg 
is similar to inorganic semiconductors (e.g. for PbTe the 
value is 0.7 x 10“'*/K at 300 KI^ and positive over the 
full temperature range. 

The temperature dependence of the phonon modes can 
be described by the Griineisen parameter, which has 
an average of around 1.6 (see Figure |^, slightly be¬ 
low the value of 1.7 found in Pblj!^. The imaginary 
modes at R and M remain at all temperatures, consis¬ 
tent with the cubic lattice being a dynamic average of 
a locally distorted structure; the same phenomenon is 
observed in GsSnlg.!^ The high level of anharmonicity 
associated with the soft titling modes is consistent with 





























the ‘ultra low’ (< 1 at 300 K) lattice ther¬ 

mal conductivity reported for single crystals and poly¬ 
crystalline MAPbIgfSS Hybrid halide perovskites are thus 
also promising for application in thermoelectric devices 
if thermal stability issues can be overcome!^ 

IV. CONCLUSIONS 

The vibrational frequencies of three crystallographic 
phases of the hybrid perovskite CH3NH3Pbl3 have been 
investigated. We identified three main phonon branches 
present in the three phases. Two high-frequency 
branches are associated with the vibration and bond 
stretching of the molecular cation with frequencies in 
the range 300 cm“^ to 3300 cm“^. The lowest en¬ 
ergy branch, below 150 cm“^, arises predominately from 
the inorganic cage, but with half the modes coupled to 
the motion of the molecule. The simulated Raman spec¬ 
trum is in good agreement with measurements on a sin¬ 
gle crystal of MAPbl3. Dynamic instabilities occur at 
the zone boundaries, which requires methods beyond the 
harmonic approximation, such as self-consistent phonon 
theory, for an accurate treatment. These results suggest 
that the room temperature structure of MAPbl3 is fluc- 
tional, owing to the persistent titling and distortion of 
the octahedral networks and rotations of the molecular 
cations. 

These factors may be important for developing a quan¬ 
titative understanding and model of how hybrid per¬ 
ovskite solar cells operate. Upon excitation, the rela¬ 
tive stability of free carriers and excitons depends inti¬ 
mately on the dielectric screening of the material, which 
includes vibrational and rotational components. The 
transport and recombination of photo-generated charge 


carriers will also be influenced by electron-phonon cou¬ 
pling, which can significantly reduce the effective size and 
distribution of electrons and holes within the perovskite 
layer. 


V. WEB ENHANCED 

An animation of the 36 P-point phonon modes of the 
cubic phase of MAPbl3 in gif format is available at http: 
//people.bath.ac.uk/aw558/temp/mapi_phonon.gif 
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